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Abstract 

A comprehensive analysis of the morphology of the solution space 
for a special type of quadratic quaternion equation is presented. This 
equation, which arises in a surface construction problem, incorporates 
linear terms in a quaternion variable and its conjugate with right and 
left quaternion coefficients, while the quadratic term has a quaternion 
coefficient placed between the variable and its conjugate. It is proved 
that, for generic coefficients, the equation has two, one, or no solutions, 
but in certain special instances the solution set may comprise a circle 
or a 3-sphere in the quaternion space H. The analysis yields solutions 
for each case, and intuitive interpretations of them in terms of the 
four-dimensional geometry of the quaternion space H. 
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1 Introduction 

The real quaternions, discovered by Hamilton in 1843, form the first known 
algebra that involves a non-commutative product, denoted by H. This prop¬ 
erty makes the solution of equations featuring quaternion coefficients and 
unknowns a much subtler and richer problem than in the case of real or com¬ 
plex numbers. The present study is devoted to developing a comprehensive 
solution for a novel type of quadratic equation, of the form 

xvx* + XQ + nx* = S, ( 1 ) 

in which both the coefficients P, Q, P, S and the variable X are quaternions 
{X* being the conjugate of <T). Since we are only interested in the genuinely 
quadratic case of equation (1), we assume henceforth that P 7^ 0. We further 
emphasize that only the algebra HI of real quaternions is considered herein: 
the results do not apply, for example, to the complexified algebra El (g) C. 

A complete analysis of equation (1) reveals that, in addition to cases with 
(at most) two distinct point solutions, special values of the coefficients yield 
singly-inhnite {circular) and triply-infinite {3-sphere) families of solutions. 
A proper identihcation and treatment of these degenerate cases is therefore 
an essential feature of any comprehensive solution procedure. 

The motivation for studying equation (1) arises [8] from the construction 
of a surface patch x(m, v) dehned on {u, v) G [0,1]^ with prescribed boundary 
curves, such that the v = constant isoparametric curves are all polynomial 
Pythagorean-hodograph (PH) curves [5]. A brief review of the surface con¬ 
struction problem can be found in Section 3 below. At present, we highlight 
some unusual properties of equation (1) that distinguish it from prior studies 
of quadratic (and higher-order) quaternion equations. 

The earliest investigations of quaternion equations are found in the papers 
of Niven [17, 18], Eilenberg and Niven [3], and Gordon and Motzkin [12]. A 
special case is that of equations of the form f{X) = 0, where f{X) is an 
element of the algebra of quaternion polynomials ]H[[A’]. In this algebra, a 
convention is fixed for the relative position of the quaternion coefficients and 
powers of the quaternion variable in each monomial. Thus, one speaks of a 
left/right quaternion polynomial if the coefficients all appear to the left/right 
of the corresponding powers of the quaternion unknown. 

Several studies [6, 13, 15] have been specihcally concerned with (monic) 
quadratic quaternion equations specihed by unilateral coefficients. In this 
case an essentially closed-form solution, requiring a determination of the 
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positive real root of a real cubic equation by Cardano’s method, is possible 
(including a complete enumeration of special cases). Although no closed- 
form solution for the roots of higher-order polynomials is possible, there has 
been considerable progress in elucidating their fundamental nature, and in 
developing numerical methods to compute them [1, 2, 9, 10, 11, 14, 16, 19, 
20, 21, 22, 23]. In particular, it has been shown that the set of roots of any 
polynomial f{X) G ]H[[A’] is a hnite union of singletons and 2-spheres. 

On the other hand, one may consider the class E of functions / : H — )■ H 
dehned as hnite sums of monomials of type AqXAi ■ ■ ■ A An with • • •, An G 
HI and n G N. This class is extremely large. In fact, it includes each of the 
functions that map any quaternion A = xq + Xii + X 2 j + Xsk to one of its 
four real components Xq, ..., X 3 since, by direct computation, 

xo = -(A — iAi—jAj — kAk), 

Xi = -^—A + iA i — j A j — kAk), 

X2 = ^{-A - iAi + j Aj -kAk), 

X3 = ^{—A — iAi—jA’j + kA’k). 

Consequently, the class of equations f{A) = 0 with / G E coincides with the 
class of all systems of four real polynomial equations in the four real variables 
Xo,. . . ,X3. Moreover, E includes the map A 1—)■ A*, so equation (1) also has 
the form f{A) = 0 with / G E. 

In the present paper, we show that the subclass of equations of the form 
(1) is special enough to allow a comprehensive solution. We show that gener- 
ically there are two, one, or no solutions, which may be determined geomet¬ 
rically. The remaining special cases, in which an entire circle or 3-sphere of 
solutions may occur, are also studied in detail. The method adopted herein 
is to reduce all instances of equation (1) to two special cases, which may be 
treated as systems of real quadratic equations — as one may expect from the 
preceding discussion concerning the class of functions E. 

The plan for this paper is as follows. First, some preparatory notations 
and results are presented in Section 2, and the surface construction problem 
resulting in equation (1) is briefly reviewed in Section 3. A comprehensive 
solution procedure, including treatment of all special cases, is then developed 
in Section 4, and examples are presented to illustrate the different solution 
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morphologies that may arise, depending on the nature of the coefficients. 
Finally, Section 5 assesses the principal results of the present study. 

2 Notations and preliminary results 

Before proceeding to the study of equation (1), we £x some notations and 
briefly recall the dehnition and basic properties of the algebra of quaternions 
H. It is the vector space endowed with a multiplicative operation through 
the following steps: 

• denote the standard basis of as 1, i, j, k; 

• let the (left or right) multiplication by 1 have no effect on l,i,j,k; 

• set F = = k^ = —1 and i j = —j i = k, j k = —k j = i, k i = —i k = j; 

• extend the multiplication to all quaternions X = xqI + Xii + X 2 ^ + Xsk 
(with xq, ..., X 3 G M) in a bilinear fashion. 

The resulting algebra is associative and non-commutative. It is also unitary, 
with identity element 1. In view of this, we henceforth write Xq instead of 
Xq 1, and we identify the axis M1 with the real held M. 

Moreover, H is a skew held — namely, every non-zero element admits a 
multiplicative inverse. In order to give a formula for the inverse, we introduce 
some further notations. For each quaternion X = xq + Xii + X 2 } + Xsk G H, 
we denote by 

X = xii + X2j + Xsk 

its vector part, and by xq its scalar part. The notations scal(T’) = xq and 
vect(T’) = xii + X 2 j + xak are also used. A quaternion whose scalar part 
vanishes is called a pure vector quaternion. To each quaternion X = xq + x = 
xo + xii + X 2 j + xsk, one associates the conjugate 

A* = xo — X = xo — xii — X2j — Xsk, 

such that XX* = Xg + |xp = Xg + x^ + X 2 + x| coincides with the square of 
the Euclidean norm \X\ of X, i.e.. 




AA* = A* A 
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By analogy with the complex case, the (non-negative) real number |A:’| is also 
called the modulus of X. Clearly \X\'^ is a positive real number for all non¬ 
zero A" G H, and the preceding equation implies that is 

the (left and right) inverse of A, i.e., 

V —1 V*l Vl~2 I vl~2v* 

A. — A, \A I — \A I A . 

The following properties will also prove useful in our computations. 

• If we denote by ( , ) and x, respectively, the Euclidean scalar product 
and the vector product in M^, which we may identify with vect(]H[) = 
Mi + Mj + Mk, then for all x, y G vect (H), we have 

xy = - (x,y) + X X y. 

• If we denote by ( , ) the Euclidean scalar product in M"^ = H, then for 
all A, 3^ G HI we have 


(A, 3 ^) =scal(A 3 ^*). 

Finally, we denote by § = {A : A^ = —1} = {x G vect(]H[) : |x| = 1} the 
2-sphere of pure vector quaternion units. Then the following properties hold. 

• For any hxed u G S the subalgebra of H generated by 1 and u, namely 
M + uM, is isomorphic to the complex plane and is denoted by Cu. 

• Cu = Cv if and only if u = ±v; otherwise Cu fl Cv = M. 

• Every A G El \ M (i.e., A = xq + x with x 7 ^ 0) belongs to Cx, where 

X 

X := -—-. 

|x| 

Every A G M belongs to Cu for all u G S. Therefore: 

H = IJ Cu. 

uES 


These properties will be used extensively in the proofs of our main results. 
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3 Surface construction problem 

The motivation for the study of equation (1) arises [8] in the construction of 
a surface patch x(m,u) for {u,v) G [0,1]^ with prescribed boundary curves, 
such that the v = constant isoparametric curves are Pythagorean-hodograph 
(PH) curves.^ Such a surface is obtained by integrating the expression 

x„(m,u) = A{u,v)iA*{u,v), (2) 

where A(u, v) is a bivariate tensor-product quaternion polynomial 

m n 

'4(«. v) = Y.Y.-^i . (3) 

i=0 j=0 

expressed in terms of the Bernstein basis 

= Q (1 - . k^0,...,d. 

The simplest non-trivial solutions correspond to m = n = 2. Integrating (2) 
allows the surface to be expressed in terms of Bezier control points pij [4] as 

x(u, ■y) = Pij bfiu) 6|(u). (4) 

i=0 j=0 

The points pij for z > 0 can be expressed in terms of the coefficients Aij for 
0 < f, J < 2, while poj for j = 0,..., 4 amount to free integration constants 
that specify x(0,u). Note, in particular, that pio and pa for z = 0,...,5 
depend only on the coefficients Am and Am for z = 0,..., 2, respectively. 

Now *^ 10 ) ^20 and Aq 2 i -^ 12 , -^22 may be used to fix the boundary PH 
curves x(zz, 0) and x(zz, 1) — i.e., to determine p^o and Pi^ for z = 0,..., 5 
— as Hermite interpolants [7]. The remaining coefficients ^ 01 ) *^ 11 ) ^21 must 
then be used to achieve desired positions for the three interior control points 
Psi, P 52 , P 53 that specify x(l,z;). This problem is under-determined, with 
three free parameters, but we assume that they are chosen a priori. One can 
formulate psi, P 52 , P 53 as quadratic expressions in ^oi )*^21 and known 
quantities. Two of the resulting equations are linear in ^ 01 ) *^21 and can be 
used to express these unknowns in terms of An. Finally, substituting these 
expressions into the third equation yields the form (1), where P, Q, 77, S are 
known, and we set X = ^n. Complete details may be found in [8]. 

^The PH curves have rational unit tangents, polynomial arc length functions, and many 
other attractive properties — complete details may be found in [5]. 
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4 The solution procedure 

Different approaches to equation (1) are appropriate, according to whether 
or not the coefficient V lies on the real axis M. These cases will be analyzed 
in detail in Sections 4.1 and 4.2 below. The principal results of this analysis 
are sketched in the following theorem. 

Theorem 1. For generic P, Q,77, iS G HI with V ^ 0, equation (1) has 
two, one, or no solutions. However, there are two non-generic instances of 
V, Q, TZ, S in which equation (1) admits infinitely-many solutions — namely, 
a 3-sphere of solutions or a circle of solutions. 

We now proceed to investigate individually the two cases P G M and P ^ M 
of equation (1), resulting in Theorems 2 and 3 below — which immediately 
imply Theorem 1. 


4.1 The case P G M 

We begin with the special case of equation (1) in which the coefficient P lies 
on the real axis M. Since we will have occasion to refer individually to the 
scalar and vector parts of the quaternion variable and coefficients, we recall 
from Section 2 their splitting denoted by T = xq + x, P = po + P; Q = % + 

P = Tq + r, 5 = Sq + s. 

Theorem 2. IfV = Po G M with po 7 ^ 0, equation (1) becomes 

Po\x\^ + XQ + nx* = S, (5) 


whose solutions in El are the points X = y — {Q* +P)/2po, where y satisfies 


|3^P=P, p: = 


scal(5) |Q + P 


* 12 


Po 


+ 


and 


4Po 

UQ 

Po 

In particular, we may distinguish the following cases. 


vect(T (Q - P*)) = vect ( S + 


( 6 ) 


( 7 ) 


1. When p < 0, there is no solution. 
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2. When p = 0, either X = —{Q* + lZ)/2pQ is the unique solution or there 
is no solution, depending on whether or not 


Po 

3. When p > 0, then 

(a) if Q = TZ* the set of solutions is the 3-sphere (6) or the empty 
set, depending on whether or not S G M; 

(b) if Q there are two, one, or no solutions, namely, the points 

X = y — {Q* + 7Z)/2pQ where y is a point of intersection of the 
3-sphere (6) with the affine line (7). Specifically, 



where A is given by 


1 

* 

to 

1 

vect (s+^^\ 


\ Po J 


A positive, zero, or negative A identifies cases with two, one, or 
no solutions. 


Proof. Setting y := X+{Q*+7Z)/2pQ and substituting X = y—{Q*+7Z)/2pQ 


into equation (5), we obtain 



= po|yr + J’A^ 



Q* + 7^p 
4po 


+27^Q+ |7^p 
2po 


Po\y\^ + yQ-Q*y* - 



-TZQ + Q*n* 
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where Q = \{Q — TZ*). Hence, setting 


5 


Q* + 7^|2 
4po 


TZQ - Q;n* 

2po 


we have 

Po\y? + yQ-&^y* = s. 

Writing 5 = Sq + §, eqnation (8) is eqnivalent to the system 


Po = So, yQ - = s. 


( 8 ) 

(9) 


The solntions to this system are those points that simnltaneonsly satisfy 
|3^|2 = P, with 

So So IQ* + 

P — — —-^--) 

Po Po ^Po 

and 

vect(3^(Q — 7?.*)) = 2vect(3^Q) = s = vect f *5 + ] . (10) 

V Po J 

If p < 0, then there is no solntion in H, whence case 1. If p = 0, then 0 is 
the only solntion to \y\‘^ = p. It is also a solntion of eqnation (10) if, and 
only if, s = 0. This verihes case 2. 

Consider now case 5, in which p > 0. If Q = 7^*, the solntion to (10) 
is either the entire space H or the empty set, according to whether or not 
s = 0. On the other hand, the solntions to (10) comprise an affine line if 
Q 71*. In this case, the solntions to (9) are the points 


3^ 







where ^ G M is a root of the real qnadratic eqnation 




TZQ 

Po 


p\Q-n*\^ 


A positive, zero, or negative discriminant leads to two, one, or no solntions. 

Finally, the translation X = y — {Q* + 7Z)/2pQ yields the solntions of 
eqnation (1) in terms of those of (9). □ 


We now describe some examples that serve to illnstrate all the possible 
types of solntion sets to eqnation (1) covered by Theorem 2. We begin with 
two examples that illnstrate cases 1 and 2 of Theorem 2. 
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Example 1. If "P = 1 = S and Q = TZ = 0, equation (1) becomes 

|A’|2 + 1 = 0, 

which clearly has no solution in H. 

Example 2. liV = Q = TZ = l = —S, equation (1) becomes 

lA’P + A’+ A”* + 1 = 0, 

which is equivalent to Iff + Ip = 0. Thus, df = —1 is the only solution in H. 

The following two examples are instances of case 3(a) in Theorem 2. 
Example 3. liV = Q = Tl=l and 5 = 0, equation (1) becomes 

|T’p ++ = 0, 

which reduces to |<T + Ip = 1. The set of solutions in El is therefore the 
3-sphere of radius 1 centered at —1. 

Example 4. liV=Q = TZ = l and 5 = i, equation (1) becomes 

IT'P + T^ + T’* = i. 

Since this cannot be satished by any fh, there are no solutions in El. 

Next is a family of examples corresponding to case 3(b) in Theorem 2. 
Example 5. If P = Q = 1 = — P and 5 = 1 + s, equation (1) becomes 

IT’P + T’-T’* = 1 + s. 

The solutions correspond to the intersections of the 3-sphere [ffp = 1 with 
the line X—X* = s. If s = 2 sin 9 s, there are two or one intersections, namely 
±cos(6') + sin(6') s. If |s| > 2, on the other hand, there is no intersection. 
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4.2 The case "P ^ M 

Consider now equation (1) when V is not a real number. In this case, we shall 
make use of a known result (see Section 2.1 of [7]) concerning the instance 
of (1) in which V, S are pure imaginary quaternions p, s and Q = 77 = 0. 
Recall from Section 2 that p := p/|p| and s := s/|s| are unit vectors in the 
direction of p and s, and Cp is the 2 -plane spanned by 1 and p. 

Definition 1. Let p, s he non-zero pure vector quaternions. If s ^ —p, the 
2-plane flp g is defined by 

ffp,s • Up s Cp , with Up s . y;; —. 

|p + s| 

On the other hand, if s = —p, then flp^g is defined to be the 2-plane through 
the origin that is orthogonal to p in the 3-space vect(]H[). 

Note that the 2-plane lip g bisects the angle between p and s. 

Lemma 1. Let p, s be pure vector quaternions with p 7 ^ 0. Then if s ^ 0, 
the 2-plane Ilp^g is the space of solutions V to equation 

Vp-sV = 0. (11) 

Moreover, the set of solutions to the quaternion equation 

= s ( 12 ) 

is the circle “^p^g in the 2-plane Ilp^g with center 0 and radius r = a/|s|/|p| . 
On the other hand, when s = 0, it is simply o := {0}. 

Although this result is already known [7], a proof is included below to make 
the presentation self-contained. 

Proof. Clearly, if s = 0, the unique solution is A = 0. We therefore consider 
only the case s 7 ^ 0. Equation (12) then implies that 



Therefore, any solution must have the form X = rlA with = |s|/|p| and 
\U\ = 1. Substituting into (12) then gives Uf)U* = s, whose solution set is 
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the intersection of the 3-sphere \Ll\ = 1 with the set of solutions V to (11). 
To verify Lemma 1, we show that this set of solutions is precisely Ifp^s- 
If s 7 ^ —p it is evident by inspection that V = Up g is a solution of (11). 
Moreover, if V satishes (11), the product of V with any element of Cp also 
satishes it. Thus, IIp^s = Up^g Cp belongs to the space of solutions of (11). We 
will now prove that the dimension of this space cannot exceed two. Indeed, 
an orthonormal transformation gives p = i and s = i cos (p + j sin y? for some 
(p G M and the rank of the linear map V V i — (i cos (p + j sin y?) V cannot 
be less than 2 since, by inspection, its values at V = i and V = k are linearly 
independent. Hence, the solution space of (11) coincides precisely with Hp^g. 
Otherwise, if s = —p, the solutions V = Uq + v must satisfy 

0 = V p + p V = 2no P + V p + p V = 2no P - 2 (v, p) 

i.e., we must have no = 0 and v T p. This is exactly the dehnition of Hp^g in 
the case s = —p under consideration. □ 

Remark 1. When s ^ —p, the circle ^p^g can he explicitly described in terms 
ofr= p|, Up^g, and a parameter 0 as 

^p,g = {t Up,g (cos 0 + sin (/) p) I 0 e [ 0, 27r) } . (13) 

Using again the notations X = xq + X, V = Po + Pi Q = Q'o + qj = ''"o + r, 
iS = So + s for the splittings of the variable and coefficients into scalar and 
vector parts, we are now ready to complete the study of the general equation. 

Theorem 3. Defining IZ and S = (sq, s) by 
n := {Q*V - TZV*){V - V *)-^, 

~ ^ Q*vn* + nv*Q n{v*-v)Q q*v*q nv*n* 

IP _ -p*|2 ' |p _ p*|2 |p _ p*|2 |p _ p*|2 ’ 

the solutions in H to equation (1) with V,Q,TZ,S G HI and V are given 
by 

T = z - {n - Q*){v - v*)-\ 

where Z G ^p,s satisfies 

{z,n) = (14) 

2|p| 

Specifically, the following cases may be distinguished. 
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1. If S = 0, then X = —{TZ — Q*){V — V*)~^ is the unique solution. 

2. If S is a non-zero real number, there is no solution. 

3. If S is not a real number, the following cases arise. 

(a) Ifn± g the set of solutions X is either the circle 

or the empty set, depending on whether or not so|p| = po|s|. 

(b) If IZ JL Ilp g there are two, one or no solutions, namely, the points 
X = Z — {IZ — Q*){V — V*)~^ where Z is an intersection point of 
the circle ^p^g with the coplanar affine line specified by intersecting 
Hp g with the hyperplane (14). 

Proof. Upon setting Z := X + {IZ — Q*){V — V*)~^ and substituting X = 
Z — {IZ — Q*){V — V*)~^ into equation (1), we obtain 

5 = ZVZ* - ZV{V* - V)~^{IZ* -Q) + ZQ 
-{IZ- Q*){V - V*)-^VZ* + IZZ* 

+ {iz- Q*){v - v*)-^v{v* - vyfiiz* - Q) 

-{IZ- Q*){V - V*)-^Q - IZ{V* - V)-\IZ* - Q) 

= ZVZ* + Z{V* - V)-^[-VIZ* +VQ + {V* - V)Q] 

+ [Q*V + IZ{-V + V- V*)]{V - V*)-^Z* 

+ \V- V*\-^{IZ - Q*)V{IZ* - Q) 

-\V- V*\-\IZ - Q*){V* - V)Q 
-\V- V*\-^IZ{V - V*){IZ* - Q) 

= ZVZ* + z{v* - vy\-viz* + v*Q) 

+ {Q*V - IZV*){V - V*)-^Z* 

+ \V- V*\-‘^{V - Q*){VIZ* - V*Q) 

+ \V- V*\~^IZ{-VIZ* + VQ + V*IZ* - V*Q) 

= ZVZ* + ziz* + izz* 

+ \V- V*\-‘^[-Q*{VIZ* -V*Q)+ IZ{-2V*Q + VQ + V*IZ*)], 
where IZ = {Q*V — IZV*){V — V*)~^. This gives 

ZVZ* + ziz* + izz* = S 


(15) 
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where 


S 


s + 
s + 


Q*{vn* -v*Q) + n{2v*Q -VQ- v*n*) 
\v — v*\^ 

Q*vn* + nv*Q , n{v*-v)Q q*v*q 

|-p _ p*|2 ' |-p _ -p*|2 |-p _ p*|2 


nv*Ti* 

|-p _ p*|2' 


Equation (15) is equivalent to the system 


Po + 2iZ* + iZZ* = So, Z-pZ* = s. 


(16) 


If s = 0, the only solution to the second of these equations is ^ = 0. This 
also satishes the hrst equation if and only if sq = 0. We have thus established 
cases 1 and 2. For case 5, with s 7 ^ 0, equations (16) are equivalent to 

ZiZ* +1ZZ* = ^°IPl ZpZ* = s. (17) 

IpI 

By Lemma 1, the solutions of the latter equation comprise the circle g of 
radius -^/Isl/lpl and center 0 in the 2-plane Llp^g. The solutions to the system 
(17) are determined by intersecting this circle with the set of solutions to 


(Z,7Z} 


go|p| -Po|§ 

2 IpI 


For any Zq G this equation is equivalent to 


{z-Zo,n) = 0 . 

Hence, the set is an affine space orthogonal to 77. Case 3(a) is verihed 
by observing that, since ^p s is a circle centered at 0 in the plane Hp g, the 
following are equivalent: 

• ‘ifp^s is the set of solutions to system (17); 

• Hp^g C Jf’; 

• 0 G and the equality {Z,1i) = 0 holds for all Z G Hp g; 

• So IpI = Po |s| and 77 T Hp g. 










15 


Finally, consider case 3(b) with TZ )L Ilp^g (whence TZ 7 ^ 0). Then Jif is a 
hyperplane that, by Lemma 1, intersects IIp^s D g in the affine line dehned 
by the eqnations 


Z p — s Z = 0, 


ziz* + izz* 


soIpI 

— 

Po|s| 


IpI 



Finally, the translation X = Z — (TZ — Q*)(P — V*) ^ gives the solntions of 
qnation (1) in terms of those of (15). □ 


Cases 3(a) and 3(h) of Theorem 3 are characterized by whether or not TZ 
is orthogonal to the plane Ifp g. In an algorithm, this conld be determined by 
checking to see if the scalar prodnct of TZ with any two linearly independent 
qnaternions in lip § vanishes. Note also that the compntation of the points Z 
in case 3(b) can be performed by: (a) writing a parameterization ^ ^ AT,+B 
of the affine line specihed by intersecting IIp^s with the hyperplane (14); and 
(b) solving the real qnadratic eqnation + B\^ = |sp/|pp. 

We now illustrate cases 1. and 2. of Theorem 3 by the following examples. 
Example 6. If P = i, Q = TZ = 1, and 5 = 0, equation (1) becomes 


XiX* + X + = 0, 


which has the unique solution df = 0 in H. 

Example 7. If P = i and Q = TZ = S = 1, equation (1) becomes 

XiX* + X + X* = 1, 


which has no solution in H. 


The following two instances exemplify case 3(a) in Theorem 3. 
Example 8. If P = i, Q = TZ = 1, and S = — i, equation (1) becomes 

XiX* + X + X* = -i. 

in this case, the set of solutions in El is the circle 


^i,-i = {x 2 j + x^k \xl + xl = l}. 
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Example 9. If "P = i, Q = TZ = 1, and iS = 1 — i, equation (1) becomes 

XiX* + x + a:* = l- i, 

which is impossible to satisfy, since the (purely imaginary) circle does 
not intersect the hyperplane X + X* = 1. 

Case 3(b) of Theorem 3 is illustrated by the following family of examples. 
Example 10. If P = i, Q = TZ = 1, and iS = Sq + i, equation (1) becomes 

XiX* + X + X* = so + l 

The solutions are the intersections of the circle i = {cos0 + sin0i : 0 G M} 
with the hyperplane X + X* = sq. If Sq = 2cos0o there are two or one 
solutions in H, namely cos0o ± sin0oi. If |so| > 2, on the other hand, there 
is no solution in H. 

Further examples of numerical solutions to equation (1), computed in the 
context of the surface construction problem (Section 3), are presented in [8]. 

5 Closure 

Although the solution of equations in the space of quaternions H has recently 
attracted considerable attention, most studies have been restricted to the case 
of unilateral coefficients. In the present study, we have considered a special 
quadratic quaternion equation in the quaternion variable and its conjugate, 
with mixed coefficients. This equation, arising from a surface construction 
problem [8], was shown to admit a complete characterization of its solutions, 
for all possible instances of the coefficients. In addition to point solutions, 
circles or 3-spheres of solutions are observed — as distinct from the case 
of unilateral coefficients, which admits [10, 19, 23] only point solutions and 
2-spheres of solutions. We have thus determined a significant class of low- 
degree quaternion equations (disjoint from that of polynomial equations with 
unilateral coefficients), for which a comprehensive solution can be achieved. 
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